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Abstract. The equilibrium equations for an isotropic Kirchhoff rod are known to 
be completely integrable. It is also known that neither the effects of extensibility and 
shearability nor the effects of a uniform magnetic field individually break integrability. 
Here we show, by means of a Melnikov-type analysis, that, when combined, these 
effects do break integrability giving rise to spatially chaotic configurations of the 
rod. A previous analysis of the problem suffered from the presence of an Euler- 
anglc singularity. Our analysis provides an example of how in a system with such 
a singularity a Melnikov-type technique can be applied by introducing an artificial 
unfolding parameter. This technique can be applied to more general problems. 
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1. Introduction 

The geometrically exact static equilibrium equations for a uniform symmetric (i.e., 
transversely isotropic) elastic rod are well known to be completely integrable [9]. In fact, 
there is a close relationship between these equations and those describing the dynamics 
of spinning tops [12]. It is also known that some perturbations of the rod equations are 
integrable, but that others are not. For instance, anisotropy of the cross-section [15] 
and intrinsic curvature [2] destroy integrability, as does the effect of gravity [6], but 
extensibility and shearability [18] do not, nor does the effect of an external force due to 
a uniform magnetic field [16]. 

In this paper we show by means of a Melnikov-type analysis that although 
extensibility and magnetic field individually do not destroy integrability, their combined 
effect does. The results may be relevant for the study of (localised) spatial configurations 
of electro dynamic space tethers, i.e., conducting cables that exploit the earth's magnetic 
field to generate thrust and drag (Lorentz) forces for manoeuvring [1,19]. 

The problem was studied by one of us in a previous paper [17], where a two- 
degrees-of-freedom system was obtained by reducing the rod equations in Lie-Poisson 
form to a canonical Hamiltonian system in terms of Euler angles and their canonical 
momenta. Numerical evidence of spatially complex rod configurations was given in the 
form of chaotic Poincare sections. A Melnikov analysis was also attempted but this 
was inconclusive as the version of Melnikov theory used, the Hamiltonian extension to 
two-degrees-of-freedom systems by Holmes & Marsden [8], was not valid due to the 
polar singularity associated with Euler angles. This singularity prevents the partial 
differentiability of the Hamiltonian with respect to the action variable that is required 
in the symmetry reduction from the two-degrees-of-freedom Hamiltonian system to a 
periodically perturbed planar system where the classical Melnikov result [7, 14] can be 
applied. 

One way around the singularity would be to use a different set of coordinates; for 
instance the Deprit-Andoyer coordinates used in several of the other studies [11,15]. Here 
we use a different method that removes the singularity from the unperturbed equilibrium 
with homoclinic orbits by introducing an artificial unfolding that corresponds to a change 
of coordinates. By a careful scaling of the equations, the magnetic rod is regarded as 
the unperturbed problem and extensibility/shearability is viewed as the perturbation 
(the opposite to the scaling used in [17]). We then apply a version of Melnikov's method 
due to Lerman and Umanskii [10] and Wiggins [22], briefly reviewed in the Appendix 
for the reader's convenience, to show that there exist transverse homoclinic orbits to a 
saddle-focus. The existence of such orbits implies via a classical result of Devaney [3] 
that the equilibrium equations for an extensible conducting rod are nonintegrable. 

Apart from its physical interest our result is also of interest because it gives an 
example of how a Melnikov analysis can be applied in the presence of an Euler-angle 
singularity. Such singularities arise naturally in many mechanical problems, especially 
in canonically reduced Lie-Poisson systems, and often coincide with the saddle-type 
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solution to which homochnic orbits are sought. In our Melnikov apphcation the 
singularity is circumvented by introducing an artificial unfolding that does not affect 
the physics of the problem. 

This paper is organised as follows. In Section 2 we briefly review the canonically 
reduced Hamiltonian derived in [16] and introduce the nondimensionalisation that will 
be used in subsequent sections. In Section 3 we introduce our scaling of the problem and 
identify the unperturbed system and its homoclinic orbit as well as the perturbation. 
The Melnikov analysis is then carried out in Section 4 and numerical illustration of the 
analytical result is given in Section 5. Finally, a summary and some comments are given 
in Section 6. 



2. Hamiltonian for an extensible rod in a uniform magnetic field 

The equilibrium equations for an extensible and shearable rod in a uniform magnetic 
fleld were derived in [17]. The equations were obtained by reduction of a nine- 
dimensional system in Lie-Poisson form with three Casimirs to a six-dimensional 
canonical Hamiltonian system in terms of three Euler angles (6*, ip, </>) and their canonical 
momenta {p0,p^,p^), similar to the reduction of the isotropic elastic rod [20] or the 
symmetric heavy top [5]. In the case of a transversely isotropic rod (i.e., one having 
equal bending and shear stiffnesses in the two principal directions of the rod's cross- 
section) the Hamiltonian is given by 

H{9,^P,<P,pe,p,,p,) = + ^ {^^'^£f^) +^,00.9 (1-1) cos^ + 



C2 (1 - j) cos^ + 1) sin e cos -^2(71 - C| - 2\p^ 

''''' ^ ^^^1 " " '^^^^^ ~ 1^^- 

Here B is the bending stiffness, J the shear stiffness, K the axial stiffness and A = IBm 
the magnetic parameter with I the current in the rod and the strength of the 
magnetic fleld, assumed uniform. Ci and C2 are constant values of two of the Casimir 
flrst integrals of the system. Ci is related to the magnitude of force in the system, while 
C2 is the force in the direction of the magnetic fleld. These constants are typically flxed 
by the boundary conditions. 

Note that is a cyclic variable and hence (the twisting moment in the rod) a 
constant. If, in addition, either the rod is inextensible and unshearable (1/J=1/-R' = 0) 
or the magnetic fleld (or electric current) is zero (A = 0) then Hamilton's equations 
corresponding to (1) have the flrst integral 

F = C2P^ + \Bcos9 - ^2C, - CI - 2\p^ (^pgsini: - cos^ ( ^'^ ~ f^^g ^ ) ) ' (2) 
rendering the system completely integrable [16]. 
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Pe - P^ij BC2 /II 



We use the constants C2 and p^ to introduce dimensionless quantities by setting 

1 ^ 

_ B\2Ci - CI) _ XB^ _ XB 

^ pi ' p'y Jp^' 

so that the dimensionless Hamiltonian H = HB/p"^ and integral F = F/{C2P^) become 

1 1 / ^ CO ^ \ ^ ^ 
H (9, ijj, pg, p^) = + 2 ( sinT'^ ) + " ^ + | ^ ^ 

+ (1 + 7 cos 9) sin 9 cos ip^J /i — 2vp^ 

+ — sin^ 9 cos^ ^{1^- ^up^) - ep^ (3) 

and 

17(0 I - - \ - /3 v//i-2t/p^ [ . /l-p^^cos^\\ 

F{9,ij,pg, p^) =p^ + -cos9 - pg sm ip - cosij [ — . (4) 

a a \ \ sm9 J J 

We shall focus on the case a = 0(1) and treat fi, u, e and 7 as small parameters in 
a scaling to be made precise in the next section (corresponding to a rod in a weak 
magnetic field with extensibility/shearability as the perturbation). 

In the absence of a magnetic field the parameter /i equals zero, provided one 
(re)defines Euler angles relative to the constant force vector in the rod, so that the 
polar singularity ^ = corresponds to the line of force. (In the reduction mentioned 
above the Euler angles were quite naturally defined relative to the direction of the 
magnetic field, which, unlike the direction of force, is constant.) Both Hamiltonian and 
first integral then reduce to the familiar expressions, the latter being equal to p^, the 
moment about the force vector. A nonzero /i thus corresponds to a different choice of 
coordinates taken inclined to the force vector; it has no physical meaning. It will be 
used in the first perturbation in the next section, however, to defiect the fixed point 
solution of the unperturbed system away from the polar singularity so that a Melnikov 
analysis can be applied. A similar inclination is used in [21] to unfold the fixed point. 

3. Scaling of the equilibrium equations 

Dropping overbars, we write the canonical equilibrium equations as 



p^ — cos9 1/(1 + 7 cos ^) sin 6^ COS'?/' ■yu , ^ 



9 = Pg, ijj = 2"- = sin 9 cos tp — e 

sin 9 y/fi — 2up^ a 

{p^-COs9){l-p^COs9) . nf, , n^ 

pg = o h asm6'(l + 7COS6') 

sin 9 

7 



(5) 



— (cos 9 + 7(cos^ 9 — sin^ 9)^ cos?/'-\//i — 2i/p^ sin 6^ cos 9 cos^ "ip (l^ — 2i/p^) 

pjf, = {1 + ycos9) sin 9 simply jj, — 2vp^ H — sin^ 9 cosipsinil) (/i — 2vp^). 

If e,7 = (i.e., in the inextensible/unshearable limit), these equations have the first 
integral F given in (4). 
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3.1. The case e = 7 = i/ = /i = 

Let e = 7 = z/ = /i = 0. Then the Hamiltonian in (3) becomes 

1 2 1 fp4, — COS . . 

H{0,ij,pe,p^) = -Po + -\ — — ^ — 1 +acos^. (6) 

This is the standard reduced Hamiltonian for the isotropic Kirchhoff rod [20]. We easily 
see that the (^,pe)-dynamics is independent of ip and therefore p^ is conserved. A 
Hamiltonian similar to (6) was discussed in Section 3 of [8]. In particular, it was shown 
that in the (6*, j90)-dynamics there exists a hyperbolic saddle point at the origin and 
a homoclinic orbit connecting it to itself at a single special value of p^. The saddle 
point corresponds to a periodic orbit in the full four-dimensional phase space. For (6) 
we have a similar situation. If (and only if) p^ = 1 and a > | then in the {0,pg)- 
dynamics there exists a hyperbolic saddle at the origin and a pair of homoclinic orbits 
(^,p,) = (±^^(t),±p^(t)), where 



^nt) = ±arccos(^l-^^sech^(^^^^ tjj, p^,{t) = e\t). (7) 

Along the homoclinic orbits we also estimate the variation of ip from to t as ip^(t), 
where 

ip^{t) = [ = - + arctan {V4:a - 1 tanh (|A/4a - 1 ^)) ■ (8) 

Jq 1 ~l~ COS 9 (tj 2 

3.2. The case e = 'j = u = and < ^ 1 

Let < /i ^ 1 while e = 7 = z/ = 0. Following the application of Melnikov theory 
to the nearly-symmetric heavy top in [8] one would seem to prove that there exists a 
hyperbolic periodic orbit near {9,pe,p^) = (0,0, 1) which gets perturbed to a periodic 
orbit with transverse intersections of its invariant manifolds. However, this is not the 
case. The problem is the polar singularity at ^ = as a result of which the homoclinic 
orbit only exists for a single value of the action integral, F = p^ = 1. The Melnikov 
result in [8] requires the homoclinic orbit to exist for an open interval of F values 
so that the partial derivative dH/dF exists, as required to reduce the two-degrees- 
of-freedom Hamiltonian system to a periodically perturbed planar system where the 
classical Melnikov theorem [7, 14] can be applied. 

A direct analysis for (5) with e = 7 = z/ = shows that there is not such a periodic 
orbit but instead two hyperbolic equilibria at 

{9, i>, pe, p^) = I arctan —,0,0, ^ ) (9) 



and 



{9,i\),pe,p^) = I -arctan^,7r,0, — p=^= I (10) 
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near {0,pQ,p^) = (0, 0, 1). We compute the Jacobian matrices for the right hand side of 
(5) with e = 7 = z/ = at (9) and (10), respectively, as 







'q;2 + /i 



/i 



-1 + 












and 



-1 + V^^^ 










v 

both of which have eigenvalues 



1 





1 












\ 

' + /i 



Hence, the equihbria (9) and (10) are hyperbolic saddles if a > |. We easily see 
that they converge to {9,ip,pg,p^) = (0,0,0,1) and (0,7r,0,l) as /i — )■ 0. Since the 
Hamiltonian system (5) is completely integrable, their stable and unstable manifolds 
cannot split. Applying a standard asymptotic analysis which is rigorously based on 
Gronwall's inequality (see, e.g., [13]), we show that they have one-parameter families of 
homoclinic orbits given by 

{9, ij,pe,p^) = {±9\t),ij\t) + i^o, ±p'^{t), 1) + 0{^) 

for t G (-To, To) with To = 0(1), where < V^o < Svr. 
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3.3. The case £ = 7 = and < z/ < < 1 

Let < u fi while e = •y = 0. By persistence of hyperbolic equilibria (see 
Proposition A.l), we see that there exist two hyperbolic equilibria in 0{v/^/Ji)- 
neighbourhoods of (9) and (10). Since the Hamiltonian system (5) is still completely 
integrable, they have one-parameter families of homoclinic orbits in Oivj^JJi)- 
neighbourhoods of (11) for t e (—T,T) with T < Tq and T = 0(1). 



4. Melnikov-type analysis 

We now let 7 = £7, where 7 = 0(1), and assume that < e ^Jv <^ ^JJi <^ 1. We 
can then apply the Melnikov-type technique of [10,22], briefly reviewed in Appendix A, 
with 

1 2 ^ ( Vi> ~ cos Q 



^0 = -^Pe + 7: — n — + a cos 6* + sin 6* cos V^a/zU - 2vp^ 
1 Z \ sm.6 



and 

ifi = — + -0:7 cos^ 6* + 7 sin 6' cos 6 cos'tp^J ii — 2z/p^ + ^ sin^ 6' cos^ ^ (/^ — 2i^p^) 

Note that hypotheses (HI) and (H2) are satisfled with K = F and k = t/^o. 

As stated in Proposition A.l, the two hyperbolic equilibria near (9) and (10) still 
persist but their stable and unstable manifolds may split. Taking T, Tq sufficiently large, 
we compute the Melnikov function (A. 4) as 



( dF dF 



dF 



+7- — sin6'(t) cos 6'(t) smi/j{t)J jj, - 2up^,{t) dt + 0(/i) 

= - r Jii-2up^{t) (peit) cosm + sin v;(t)') dt + 0(/.) 

= ±^i^^?^sin^o + 0(;.), (12) 
a 

where the above integrands are evaluated on the homoclinic orbits for /i, z/ 7^ and 

sin^^'(t) 



Ma) 



00 



■cosV^^(t) -p^(t) sinV^^(t) 



dt. (13) 



1 + cos6"^(t) 

Here we have used the fact that 0^{t) is an even function of t, p^(t),ip^{t) are odd 
functions and the integrand of (12) tends to zero exponentially as t — ±00. Note that 
M{iPq) is independent of 7 up to 0(^/1). 

Surprisingly, we can compute the integral (13) analytically as follows. Let 6 = 

a — i. By (7) and (8), we have 

sin^^(t) 2Ssech. St v , , 4 sech 5t tanh 5t 



l + cos^i>(t) ^/A5^tanh.^5t + l y^A6H^^^^?StTl'' 
^^(t) = it + arctan(2(5tanh(5t). 

Using the above relations and noting that 
cos (arctan(25 tanh 6t)) 



y/4:6^ tanh^^t + 1 

, e \\ 25 tanh 5t 

sm (arctan(2o tanhotj) = — ^^^^^^^^^= 

V4(52 tanh^^t + 1 




Figure 1. Plot of the function A (a). 



we can rewrite the integral (13) as 



/oo / TT \ 

sech 6t cos dt = 2n sech ( — ) 



for a > |. The function A(a) is plotted in figure 1. 

Thus, if e 7^ the Melnikov function M{iPq) has a simple zero near i/^q = 0, vr 
and by Theorem A. 2 there exist transverse homoclinic orbits to equilibria near (9) and 
(10). This result holds whether or not 7 = 0. Moreover, it follows from Devaney's 
theorem [3] that the system is chaotic in the neighbourhood of the homoclinic orbit in 
the sense that any local transverse section to the homoclinic orbit contains a compact, 
invariant, hyperbolic set on which the Poincare map is topologically conjugate to a 
Bernoulli shift of finite type. This in turn implies that the Hamiltonian H has no analytic 
conserved quantities independent of H itself, i.e., the corresponding Hamiltonian system 
is nonintegrable [15]. 

5. Numerical analysis 

Here we give numerical evidence of the above theoretical results. We first describe our 
approach to compute the stable and unstable manifolds. 
We write (5) as 

i = JDHiO, (15) 
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where ^ = {6,ip,p0,p^) and J is the 4x4 symplectic matrix, 

( idA 
[-id, ) ' 

with id2 the 2x2 identity matrix. Let denote the equihbria near (9) and (10). We 
compute W^{C,o) as follows. 

We approximate W^{^o) by the unstable subspace of the linearised system 

f] = JB^H{^o)v (16) 

near = C^q. Under this approximation, we obtain an orbit ^{t) on W^{C,o) as a solution 
of (15) satisfying the boundary conditions 

L,m = o, e(Tu) = c (17) 

where Tu > is a constant, Lg is a 2 x 4 matrix consisting of row eigenvectors 
corresponding to eigenvalues with negative real parts of the Jacobian matrix JD'^H{^q), 
and ^0 is a point on W^{^o). Thus, we solve the boundary value problem (15) and (17) 
and continue the solution in to compute W^"(^o) numerically. 

Similarly, to compute W^{^o), we solve the boundary value problem (15) and 

Lum = o, a-Ts) = eo, (18) 

where Tg > is a constant, Lu is a 2 x 4 matrix consisting of row eigenvectors 
corresponding to eigenvalues with positive real parts of the Jacobian matrix JD^H{C,o), 
and is a point on W^{^o)- The intersection of W^{^q) and Vr"(^o) gives a homoclinic 
orbit to .^0- We adopt a solution of the linearised system (16) with T^^,Ts small as a 
starting solution, and use the numerical continuation tool AUT097 [4] to perform the 
above computations. In these computations we can also take Tg^u, ^(0) and ^^'^ as free 
parameters. 

Figure 2 shows a numerically computed equilibrium (6*0,^^0) near (9) in (5) for 
6 = 1/ = 0.01, 7 = and a = 0.5. Note that {'4'iPe) = (0, 0) at the equilibrium and there 
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Figure 3. Numerically computed stable and unstable manifolds of the equilibrium 
near {0,ip,pg,p^) — {9q,0,0,p^o) in (5) for /i = 0.1, e ~ ly ~ 0.01, 7 = and a = 0.5: 
(a) ^ = 0; (b) pg = 0. The stable and unstable manifolds arc plotted as solid and 
broken lines, respectively. In plate (a) "•" represents the position of the equilibrium. 




Figure 4. Numerically computed transverse homoclinic orbit in (5) for ^ = 0.1, 
e = = 0.01, 7 = and a = 0.5. 



is another equilibrium at {9, ipyPeyPi^) = {—^o, tt, 0,p^o) near (10). These equihbria exist 
only for /i > 2z/ = 0.02. 

Figure 3 shows {ip = 0} and {pg = 0} sections of numerically computed stable and 
unstable manifolds of the equilibrium {9,ip,Pe,Pi}) = {9o,0,ip,p^o) near (9) in (5) for 
H = 0.1, e = u = 0.01, 7 = and a = 0.5. From these sections we see that these 
manifolds intersect transversely in the energy level set and there exists a transverse 
homoclinic orbit, as predicted by the theory in Section 4. We note that other branches 
of these manifolds, which intersect {6 = 0}, are very difficult to compute due to the 
singularity in (5). The transverse homoclinic orbit detected in figure 3 is plotted in 
figure 4. 
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6. Conclusion 

We have shown that the reduced equihbrium equations for an extensible and shearable 
conducting rod in a uniform magnetic field possess transverse homoclinic orbits to 
saddle-type equilibria. Via a classical result in dynamical systems this implies that 
the equations are nonintegrable and give rise to chaotic dynamics. To obtain the 
result, we applied a version of Melnikov's method, taking the classical Kirchhoff rod 
as the unperturbed problem and introducing magnetic and elastic effects as successive 
perturbations. 

Our method is also of interest because of the way it deals with the Euler-angle 
singularity of the reduced equations. Such a singularity is quite common in mechanical 
problems, e.g., in rigid-body dynamics. Indeed, it was present in the nearly-symmetric 
heavy top example treated in [8]. As pointed out elsewhere [21], in this example the 
conditions required for the Melnikov method are not satisfied and the application is 
therefore invalid, although there can be no doubt that the result is true (i.e., the nearly- 
symmetric heavy top is nonintegrable). In a certain limit, namely 7 = = e = 0, 
/i ^ a = 0(1), our magnetic rod gives an example where a similar cavalier application 
of the Melnikov method would suggest that the system is not integrable, yet the presence 
of the explicit first integral F implies that the system is in fact integrable. The resolution 
is that the fixed point {6,pe) = (0, 0) of the planar unperturbed system is not perturbed 
into a periodic orbit of the /z-perturbed system. Rather, two fixed points appear 
(cf. (9) and (10)), unfolding the Euler-angle singularity. As explained in Section 2, the 
parameter fi is an artifical parameter if no magnetic field is present. Thus we see that 
in our Melnikov application the singularity is circumvented by the introduction of an 
artificial unfolding that does not affect the physics of the problem. A similar technique 
(an inclined plane) is used in [21] to give a correct proof of the nonintegrability of the 
nearly-symmetric heavy top. The method may be applicable to other systems. 
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Appendix A. Melnikov-type technique 

For the reader's convenience, we briefly review an analytical technique to detect 
the existence of transverse homoclinic orbits to hyperbolic saddles in two-degrees-of- 
freedom, nearly-integrable Hamiltonian systems. See [10,22] for details. 
Consider a two-degrees-of-freedom Hamiltonian systems of the form 



where Hq, Hi : M. are C"'^^ (^ > 2). When e = 0, equation (A.l) becomes 



X = JBHoix) + eJDHi{x), x e R^ 



(A.l) 



X = JDHoix). 



(A.2) 
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We make the following assumptions on (A. 2). 

(HI) The system (A. 2) has a first integral K{x) such that Hq{x) and K{x) are 
independent, i.e., J^Hq{x) and 'DK{x) are linearly independent, and in involution, 
i.e., their Poisson bracket is zero, 

{A', H^} = DK{x) ■ JDHoix) = 0. (A.3) 

In other words, the system (A. 2) is completely integrable. 

(H2) There is a hyperbolic saddle point at x = Xq possessing a one-parameter family 
of homoclinic orbits x^(t; k), k G J^, where C M is a nonempty open interval. 
In addition, x^(t; k) is not only in t but also in /t. 

The Hamiltonian system (A.l) also has the System III form of [22] with n = 2 and 
m = 0. 

Since equation (A. 2) is Hamiltonian and four-dimensional, the hyperbolic saddle 
X = xq has two-dimensional stable and unstable manifolds, which are denoted by 
iy''^(xo). Let T = {x = x'^it; k) 1 1 G M, K G J^} U{xo}. It follows from assumption (H2) 
that 

w\xo) n wi^xo) D r. 

By a rather standard result on persistence of hyperbolic equilibria and their stable and 
unstable manifolds (see, e.g., [13]), we immediately have the following result. 

Proposition A.l For e > sufficiently small the system (A.l) has a hyperbolic 
equilibrium at x = Xe = xo + ff{e). Moreover, the stable and unstable manifolds W^''^{xs) 
of Xe are ff{e)-close to W^'^{xo). 

Now define the Melnikov function as 

/oo 
BK{x^{t; k)) ■ JBHi{x^{t; K))dt. (A.4) 
'OO 

We have the following result, as shown in [10]. 

Theorem A. 2 Suppose that M{k) has a simple zero at k = kq, i.e., 

dM 

M(«:o)=0, — (ko)^O. 

Then for e > sufficiently small the stable and unstable manifolds W^'^{xe) intersect 
transversely and there exists a transverse homoclinic orbit x = x^{t; hq) + G{e) to the 
hyperbolic saddle x^. 
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